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Summary.  The  paper  addresses  the  problem  of  a-posteriori  estimation  of  the 
modelling  error  and  the  adaptive  modelling  of  plates.  It  proposes  a  family  of 
models  and  provides  upper  and  lower  a-posteriori  estimates  of  the  modelling 
error  in  various  norms  .  The  estimates  eu'e  locally  asymptotically  exact. 

1.  Introduction 

The  reliability  of  computational  engineering  einalysis  depends  on  the 
reliability  of  the  underlying  mathematical  formulation  (i.e.  the  mathematical 
model)  of  the  engineering  problem  under  consideration  as  well  as  on  the 
reliability  of  the  numerical  treatment  of  the  mathematical  model. 

The  selection  of  the  mathematical  formulation  is  the  most  crucial  aspect 
of  reliability  in  the  analysis  of  an  engineering  problem. 

Reality  czui  never  be  completely  identified  with  any  mathematical  formu¬ 
lation.  Nevertheless,  in  practice  reality  is  identified  with  the  exact 
solution  of  a  complex  mathematical  problem  to  which  we  shall  refer  as  “true 
problem"  throughout.  Prior  to  its  computational  solution,  however,  the  true 
problem  is  approximated  by  a  structurally  simpler,  so-called  mathematical 
model  and  then  the  computational  analysis  is  performed  on  this  selected 
mathematical  model.  This  means  that  the  reliability  of  the  computed  solution 
should  not  be  measured  against  the  mathematical  model,  but  rather  against  the 
true  problem,  i.e.  it  should  take  into  account  the  discretization  and  the 
modelling  error.  (For  a  discussion  of  mathematical  modelling  in  engineering 
we  refer  to  (11). 

The  usual  engineering  approach  is  to  test  the  reliability  of  a  mathemati¬ 
cal  model  by  a  set  of  benchmark  analyses  for  which  the  solution  of  the  true 
problem  is  known.  This  leads  often  to  the  design  of  various  models  tailored 
to  different  objectives  of  the  engineering  euialysis. 

No  single  mathematical  model  ceui  be  good  (or  optimal)  under  all 
circumstances  of  the  engineering  analysis  and  therefore  a  decision  has  to  be 
made,  upon  which  particular- model  the  computational  analysis  should  be  based. 
The  strategy  of  the  model  selection  must  be  governed  by  an  a-posteriori 
assessment  of  the  modelling  error  and  an  adaptive  selection  of  the  model 
(without  knowledge  of  the  exact  solution  of  the  true  problem).  The  derivation 
of  reliable  a-posteriori  estimators  for  the  modelling  error  is  the  topic  of 
the  present  paper.  More  specifically,  we  will  discuss  the  hierau'chical 
modeling  of  thin  structures,  such  as  plates  and  shells,  in  mechanics. 
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The  Plate  and  shell  problems  and  their  hierarchical  models. 

The  problem  of  plate-and  shell  analysis  has  been  in  the  forefront  of 
engineering  interest  for  a  long  time.  Here  the  true  problem  is  the  three- 
dimensional  problem  of  linear  elasticity  on  a  thin  domain  -  a  planair  or  curved 
surface  of  thickness  d.  Plate-and  shell  models  are  then  certain  2-dlmensional 


problems,  i.e.  the  modelling  leads  to  a  reduction  of  the  dimension  of  the 
problem  by  one. 

A  particular  model  still  widely  used  today  was  proposed  already  in  the 
first  half  of  the  nineteenth  century  by  S.  Germain  [2]  Euid  G.  Kirchhoff  (3]. 
The  derivation  of  this  and  other  models  is  mostly  beised  on  physical  considera¬ 
tions  or  mathematical  analyses  of  various  degrees  of  rigor  (see  e.g.  [4]-[19] 
and  the  references  therein)  or  on  asymptotic  analysis  of  the  three- 
dimensional,  true  problem  sls  the  thickness  d  tends  to  0  (see,  for  example 
[9],  [20]-[28]  and  the  references  therein).  Other  results  address  the 
theoretical  relation  between  the  true  problem  amd  the  models  (see  e.g.  (29]  - 
[32])  or  they  investigate  the  derivation  of  hierarchical  families  of  models  ais 
general  problem  of  dimensional  reduction  (see  e.g.  [33]  -  [36]  for  related 
ideas). 

The  solutions  of  various  proposed  models  could  lead  to  significantly 
different  results  which  explain  the  very  lau'ge  number  of  proposed  formulations 
in  the  literature  today.  For  theoretical  and  niimerical  analyses  of  various 
models  see  [37],  [38].  Let  us  consider  as  am  example  for  the  difference 
between  various  models  the  claunped  squau'e  plate  :  = 

(-0. 5, 0. 5)^ X  (-0. 005,0. 005)  of  thickness  0.01  with  uniform  load.  We  will 

7 

aissume  the  Young’s  modulus  of  elaisticity  E  =  10  amd  a  Poisson  ratio  r  =  0.3. 
Let  us  be  interested  in  the  bending  moments  Mjj(Xj,0)  and  M22[5tj,0).  Table 
1  shows  the  data  for  the  three  dimensional  ("true")  formulation,  the  Reissner- 
Mindlin  model  (RM)  and  the  model  (1,1,2)  (see  [37],  [39],  [40]  for  the 
definition  of  (1,1,2)  model). 

We  see  that  the  reliability  of  the  model  strongly  depends  on  the  aim  of 
the  auialysis.  If  we  are  Interested  in  Mjj(0.4,0),  M22(0.4,0)  then  the 
diff^i-ence  between  the  models  is  practically  negligible  but  in  the  caise 
M22(0-S,0)  the  difference  is  significamt  (.307.) 
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Table  3. 1. 


The  flioments  and 


for  vzirious  models. 


Mii(Xi,0) 

M22(Xi,0) 

3  dim 

RM 

(1.1,2) 

3  dim 

RM 

(1,1,2) 

0. 

-0.0229 

-0. 0229 

-0.0229 

-0. 0229 

-0. 0229 

-0. 0229 

0.2000 

-0.0157 

-0.0157 

-0.0157 

-0.0163 

-0.0163 

-0.0163 

0.4000 

-0.0164 

-0.0163 

-0.0164 

-0.0026 

-0.0027 

0.0026 

0. 4900 

0. 0470 

0.0164 

0. 0470 

0.0141 

0.0141 

0.0141 

0. 4990 

0.0509 

0.0509 

0.0509 

0.0168 

0.0152 

0.0170 

0.4993 

0.0510 

0.0510 

0.0510 

0.0176 

0.0153 

0.0179 

0. 4999 

0.0512 

0.0512 

0.0512 

0.0207 

0.0153 

0.0212 

0.5000 

0.0513 

0.0513 

0.0513 

0.0220 

0.0153 

0.0220 

3.  The  bovmdarv  value  problem  aind  the  hierarchy  of  the  models. 

To  present  all  the  essential  ideas  of  adaptive  modelling  in  the  simplest 

setting,  we  will  restrict  ourselves  to  the  problem  of  the  Laplace  equation  on 

a  thin,  plate-like  domain.  Nevertheless  the  approach  presented  here  can  be 

used  in  general  [41]. 
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Let  (i)cR  denote  a  bounded  domain  with  piecewise  smooth  boundary  P 

d  3 

8Uid  define  the  plate  of  thickness  d,  fl  :=  {x  =  (Xj,X2,X2) e R  | (Xj,X2) € w, 
IXgl  <d/2>.  Further,  we  introduce  the  lateral  boundary  S  and  the  faces  of 
the  plate  R^: 

S  :=  {x  e  R^l (Xj,X2) €  r, IXgl  <  d/2},  :=  {x €  R^| (x^, X2) e  u, x^  =  ±d/2}. 

We  will  be  interested  in  the  problem 


Aw  =  0 

on 

d". 

(3.1) 

o 

II 

on 

s. 

3w  _  1 
an  ■  2 

on 

2 

Here  f(Xj,X2)eL  («}  and  n  denotes  the  exterior  unit  normal  vector. 


Let  us  cast  (3.1)  into  the  weeik  form.  Denote  by 

(3.2)  H  :=  {u€H^(n^)|u  =  0  on  S} 

and  define  by 


lync  QUALrrf  inspected  3 


.  - 

j  D I  s»  on/ _ 

i  4vellatinty  CodM 
j  j  A  vail  aoidicr 

[Diet  I  Speal^ 

i 


□a 


Obviously  here  the  solution  Is  symmetric,  i.e. 

C3.6)  u(Xj,X2,X2)  =  uCXj.Xg.-Xg). 

Assuming  that  the  thickness  d  of  is  small,  we  will  approximate  the 
full,  three  dimensional  problem  (3.1)  -  (3.6)  by  a  sequence  of  two  dimensional 
problems,  the  hierarchical  family  of  plate  models.  Let  us  define  these 
models. 

Denote  by  f  =  (Wj^  £  w,  l^i^n)  a  partition  of  u  into  n  domains 

w.  with  piecewise  smooth  boundaries  du.  such  that  u.r\u  =  ^  for  i*J 
_  n  _ 
atnd  u  =  U  u. . 
i=l 

We  will  allow  the  order  of  the  models  to  be  different  in  different 
subdomains  euid  define 


<I  :=  integers. 

Denoting  by  Lj  the  Legendre  polynomial  of  degree  J  on  (-1,1),  we  introduce 
the  space 

SCf.q)  {u.Hlul  .  ^ 

J=0 

and  define  the  (?*,  q)-piate  model  as  the  boundairy  value  problem; 

Find  u(y,  q)  €  S(?’,  q)  such  that 

(3.7)  B(u(?',q),v)  =  F(v)  Vv€S(?’,q). 
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Obviously,  due  to  (3.6),  we  can  assume  a  priori  that  ^  for  odd  J. 

The  solution  u(?, q)  of  the  (y,q)-plate  model  is  the  (energy)  projection  of 
the  three  dimensional  solution  Uq  of  (3.5)  onto  S(?’, q). 

We  have 


Theorem  3.1.  Let  e 


(3.8) 


almost  everywhere. 


UQ-u(y,q).  Then 


d/2 

1  ' 


e(Xi,X2,X3)dx3  =  0 


-d/2 


Proof  For  any  U(x.  ,x,,)  €  S(y,q)  we  have  B(e,U)  =  0.  Observing  that  =  0 

1  A  on 

on  and  integrating  by  parts  we  get  (3.8). 

Now  we  state  the  goal  of  the  computation  (and  of  the  hierarchical 
modelling) : 

Given  load  data  f,  a  norm  |ol  on  H,  and  a  tolerance  t > 0,  find  a  model 
(y, q)  so  that 


(3.9) 


|uQ-u(f,q)|  <  T. 


‘  Remark  3. 1.  If  any  re2isonable  norm  |o|,  such  that  lu^l  <  w,  and  the  fixed 
thickness  d  of  the  plate  are  given,  then  such  a  model  exists;  Select  in 
(3.7)  n  =  1  and  f  =  {u},  q  =  {q>  sufficiently  large.  a 


Remark  3.2.  We  assumed  here  that  the  (f, q)  plate  model  can  be  solved 
exactly.  This  is,  of  course,  not  true  in  practice,  since  typically  (3.7)  must 
also  be  solved  numerically  by  e.g.  the  finite  element  method.  We  distinguish 
between  the  modelling  error  (3.7)  and  the  discretization  error.  Since  we  are 
here  primarily  interested  in  the  modelling  error,  we  eissume  below  that  the 
discretization  error  is  negligible,  l.e.  that  the  (3’,q)  model  is  solved  to 
sufficiently  high  accuracy.  □ 


Remark  3.3.  All  the  notions  introduced  above  apply  also  verbatim  to  the  plate 
problem  discussed  in  Section  2.  The  model  problem  (3.1)  -  (3.7)  exhibits  the 
bsLSic  characteristics  of  the  plate-  and  shell  problems.  □ 
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Remark  3.4.  The  selection  S{.T,q)  in  (3.7)  is  not  arbitrary.  It  vras  shown  in 
[33],  [34]  that  this  selection  has,  for  fixed  two  essential  prop)erties, 
namely: 

a)  the  selection  of  polynomials  in  is  the  only  choice  which  leads  to  the 
optimal  rate  of  convergence  as  d — >0  (although  for  fixed  d>0  this 
choice  is  not  optimal  [34]), 

b)  as  — >00,  1  S  i  5  n,  u(y,q) — >Uq  and  a-priori  error  estimates  are 

given  in  [33],  [34]. 

For  laminated  beams  and  plates,  an  alternative  subspace  S(f, q)  with 
similEU*  properties  is  known  [17],  [19].  a 

Remark  3.5. 

The  adaptive  selection  of  the  model  has  now  two  essential  parts: 

i)  Design  of  the  local  error  Indicator,  which  leads  for  the  norm  under 
considerations,  to  the  error  estimator, 

ii)  Design  of  the  adaptive  procedures. 

We  will  address  these  points  in  the  next  sections.  □ 


4)  Some  abstract  results 

Let  Hj,  Hg  be  two  reflexive  Banach  spaces  furnished  with  the  norms  ll'Hj 
and  Il•ll2»  respectively.  Further,  let  B(u, v)  be  a  biline^u'  form  defined  on 
Hj X  Hg.  We  will  call  the  bilinear  form  (C, y)-reguJar  if  there  exist 
constants  0  <  C,  7  <  oa  so  that 

(4.1)  |B(u,v)|  iCIlulljlIvlIg, 

(4.2)  inf  sup  |B(u, v)|  i  y, 

l|u||j=l  ||v||2=l 

(4.3)  for  any  v*0,  v€H„,  sup  |B(u,v)|  >  0. 

l|u||j=l 


Bilinear  forms  satisfying  (4.1) 
(see  [42]). 

a)  Let  fedig)'  (i.e.  f  is  a 
there  exists  exactly  one  ueH^  such 

B(u.v)  = 

b)  If 


-  (4.3)  have  the  following  properties 

bounded,  line^u'  functional  on  H2),  then 
that 

f(v),  VveHg. 
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(4.4) 


then 


sup  |B(u, v) I  ^  A, 

I|vll2=l 


ilullj 


< 


A 

7‘ 


Let  us  consider  now  some  special  cases  which  will  be  important  later. 
Theorem  4. 1.  Let 


ue  H 


-d/2 

J 


(4.5)  =  I 

-d/2 

where  H  is  as  in  (3.2).  Let  further 


X2)dx2  =  0  almost  everywhere 


(4.6) 

llullj  = 

"  M 

r 

1 7u 1 ^dx 

■ 

■ 

ll/2 


=  u 


Then  the  bilinear  form 
(3.3)  B( 


u.v)  :»  I 


Vu  •  7v  dx 


Si 


is  (C, y)-regular  with  C  =  y  =  1. 


Theorem  4.2.  Let  9(Xj,X2)>0,  (Xj^,X2)eu  have  bounded  first  derivatives 
aind  define 

(4.7) 


^  |ll"|  /  /ix  X  ). 

(x,.x.,)€«  l‘^''il  ^  ^  ^ 


r  2 
1=1,2 


Further  let  =  H2  where  Hj  and  H2  are  the  same  spaces  as 

but,  instead  of  |l«||j  aind  Ihllg* 
respectively,  which  au'e  given  by 


before,  but,  instead  of  |1«|L  aind  ll«IU,  furnished  with  the  norms  ||o||  and 

1  ^  i,<p 


"  ”2.9 

(4.8) 

(4.9) 


J  ^ 


^IVul^dx, 


Q 


iivii: 


■  I » 

J  ^ 


■^iVvI^dx  =  ||v||^ 


l.VP 


-1' 


a 
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Then  the  bilinear  form  (3.2)  is  (C, y)-regular  with  C  =  1  and, 


(4. 10) 


riro  - 


■1/2 


Proof.  By  the  Schwarz  inequality  we  see  immediately  that  (4.1)  holds  with 
C  »  1. 

Now,  let  ucH!?,  Hull,  =1  and  define  v  =  ^^u.  Then  veH??  auid 

1  1.  <p  ^ 

2 

dv  2du  ^  dip  .  ,  „ 

dx^  “  ^  dx^  “dx^'  ^ 

dv  2  du 

Hence 


(4. 11) 


B(u.v)  =  iiuii; 

1.  V 


r  1 

r  2 

2 

f 

d/ 

J 

^  i=l 

dx 

J  1 

1 

dx. 


Q 


A/2 


Realizing  that  J  u(Xj, Xg, X2)dX2  =  0  almost  everywhere,  we  have 


-d/2 
-d/2 


(4. 12) 


d/2  2 


dx- 


-d/2 


-d/2 


Hence,  with  (4,7)  and  i  =  1,2  we  have  for  any  e>0 

2 


I  I  ‘**1  I  (I^J 

n”  o'* 

c  %  I  (%]  5  \  ''') 


O'" 


Selecting  c  »  we  get 


r  1 

r  ^ 

9  1 

1 

Y*  du 

Z-  3x7' 

^  i=l 

d/ 

dx^ 

J 

dx 

id  1 

n 


and 
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|B(u.v)|  i  Hull 


i.<pL  "  J 


similarly 


Hvllp  ^  =  llull^  ^  + 

2,<p  l.fp 


•1=1  J 


Hence 


1  ir  ^ 


-1/2 


Since  B(u, v)  Is  symmetric,  (4.3)  holds  euid  the  theorem  Is  proven.  n 

We  can  select  In  particular  any  ^(x^.Xg)  so  that  2Q  i  ^  •  a,  0<a<l, 
and  get 

1/2  -1/2 

7  S  (l-a)^^‘'{l  -••ad  +«))  c 


dtp  dtf) 

since  =  2<p  vfe  have 

dx^  ^  dx^ 


(4. 13) 


Q  :=  2  max 

(Xj,X2)€W 

1=1,2 


d(p 


ax. 


I«>l- 


Of  particular  Interest  In  sections  5  and  6  will  be  the  weight  function 


(4. 14) 


=  e 


-B(  IXj-XjdIXg-Xgl ) 


where  Q  =  2B.  Here  we  can  select  40  =  3a,  0<a<  1  and  obtain  that 

d 


B(*,*)  Is  (C, y)-regulzu’  on  H^  x  H^- 


We  note  that  for  0  <  <  ^  <  <2  ®  simple  au'gument 


y  2  — 


Let  us  now  consider  the  bilinear  form 


(4. 15) 


«  • 

defined  on  x  Hg  where 


Bj(u,v)  =  I 


uAv  dx 


n 
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H. 


-d/2 

ueL2(n)l  u(Xj.X2.X2)dx2  =  0 


almost  everywhere 


-d/2 


fvirnished  with  the  norm 


|u|j  * 


^l/2 


lul^dx 


8md  where 


(a'^l  I  Il4v|i^,„d,  <  -.  *  0  on  R, 


4 


• 

-  • 

|v|2  = 

1 Av 1 ^dx 

B  * 

.d 

where  H  is  defined  in  (4.5),  fiirnished  with  the  norm 

,1/2 

0-' 

We  note  that  locally  and  hence  meikes  sense.  It  also  follows 

^  *  an 

easily  that  1*12  Is  a  norm  on 

It  is  easy  to  see  that  (4.1)  holds  with  C  =  1.  We  will  now  estimate  7. 
For  given  ueH^  let  us  define  s  to  be  solution  of  the  equation 

(4.16)  As  =  u. 


(4. 17) 

(4.18) 


s  =  0  on  S, 


3H*°  "f 


Because  u€L2(R)f  s  obviously  exists  and  is  uniquely  determined. 

.d/2 


Define 


=  s-^  s(Xj,X2,X2)dx, 


(4. 19) 

a 

-d/2 

then  also  (4.17)  and  (4.18)  holds  for  z. 
Further,  by  using  (4.18) 

.d/2 


Az  =  As  -  ^ 
o 


(As)(Xj,X2,X2)dX2, 


-d/2 


hence 
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.d/2 

=  u  -  i  uCxj.x^, 


X3)dx3  =  u. 


Therefore 


}j(u.2)  =  I 


uAzdx  *  lAzlIlg^lllulllj 


and  (4.2)  follows  with  7=1. 

Hence  we  have  proved 

Theorem  4.2.  The  bilinear  form  B^tu.v)  defined  by  (4.14)  on  H*  x  H* 
satisfies  (4.1)  and  (4.2)  with  C  =  y  =  1. 


5.  A  posteriori  estimates  of  the  modelling  error. 

In  this  section  we  assume  that  the  (?*,q)-plate  model  is  given  and  its 
exact  solution  u(?,q)  defined  in  (3.7)  is  known.  We  will  be  interested  in 
estimating  the  modelling  error 

luQ-ulf.q)! 

for  a  particular  norm  \o\  which  is  L^-based.  With  the  norm  |o|  we 
associate  the  error  indicator  function  so  that,  for  1  <  p  <  », 


(5.1) 


®(u(?’,q))  =  |7)(x,,x  )l^dx,dx. 

1  &  1  < 

•  ^ 


»  |u3-u(3>.q)|. 


The  qxiantlty  8  will  be  called  error-estimator.  The  effectivity  index 
corresponding  to  8  is  defined  by 


(5.2) 


a  ._  S(u(f.q)) 

®  •  ~  rr: - ../a  _ 


IUq  -  u(?,i 


We  say  that  g  is  an  upper  (lower)  estimator,  if  0>1(0<1), 
respectively,  and  call  the  estimator  S  proper  with  respect  to  a 

class  T  of  data  f,  if 


(5.3) 


holds  for  all  data  f € T. 


O<KjS0Sic2<® 
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The  estimator  S  Is  asymptotically  respectively  spectrally  exact  on  T, 
if 

(5.4)  @ — >1  as  d — >0*  respectively  as  q — 


for  Uq  belonging  uniformly  in  d  to  certain  classes  of  data. 

Vfe  shall  say  that  the  Indicator  function  ijCXj.Xg)  is  locally 
asymptotically  exact  (resp.  locally  proper),  if  (5.4)  (resp.  (5.3))  holds  for 
the  norm  |o|  =  ||o||  where  the  weight  is  given  by 

l,<p  a  a  ^ 


(5.5) 


^(x^.Xg) 


:=  exp 


0, 

Xi-Xjl 


]■ 


p  e  (0, 1 ) ,  0  <  a 


and  ^Xj.Xgjew  arbitrary. 

We  will  auidress  here  for  simplicity  only  the  case  n=l.  p=2,  T  =  {w> 
and  q  “  {q},  qiO,  and  denote  the  corresponding  subspace  by  S(q). 

Throughout  this  section  we  will  denote  by  ||o|l^  ^  the  norm  defined  in 
(4.8)  auid  the  norms  of  f  are  always  understood  8is  integrals  over  w  only. 
Let  us  first  analyze  the  modelling  error  in  this  norm.  We  begin  by  zmalyzing 
the  best  possible  asymptotic  rate  of  convergence  of  the  (P,q)-model. 


Theorem  5. 1  Assume  that  qSO  is  an  integer  euid 

(5.6)  €  (fi^  nH^)(w)  for  1  i  J  S  q. 

Then  there  exists  u(2q)eS(2q),  u^  -  u(2q)  €  (see  (4.5))  such  that 

(5.7) 

Proof:  We  only  sketch  it  since  it  follows  closely  that  of  [33,  Theorem  3.1]. 
First  observe  that  for  T ^ Tq  as  in  (4.10) 

|B(u--u(2q),v) I 

(5.8)  y  l|u  -u(2q)||  5  sup  - ^ -  . 

llvll  ^  "''"2.9 

2.9 

Now,  in  [33]  it  was  shown  that  there  exists  a  sequence  of  polynomials 
^2^  degree  2J,  J  *  0,1,2,...  which  is  independent  of  d  such  that 
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G(2q) 


J=0 


satisfies,  for  every  veH, 


B(Uq  -  u(2q), v)  =  d 


2q+l 


J  J  «*3  ‘ 


-d/2  w 


where  ^2q^^l’^2^  “  ^2q+2(7dr]^^^^^^^l‘’^2^‘  (5.7)  is  a  simpl 


of  Schweu'tz’s  inequality  and  of  (5.8). 


e  consequence 
□ 


2 

Remark  5.  1  If  q  =  0  then  feL  (<*»)  is  sufficient  for  (5.7)  to  hold.  a 

Remeu^k  5.2  By  definition  (3.7)  of  the  (f.q)  model  and  Theorem  3.1  the 
approximation  is  quasioptimal,  i.e. 


(S.9)  |u  -  u(r,q)|  S  C(,  )  Irf  ° 

X€S(T,q) 

To  derive  the  error  indicator  function  tilXj.Xg),  we  note  that  using  Theorem 
3.1  the  modelling  error  eg^  :=  UQ“u(?*,q)  belongs  to  defined  in  (4.5) 
and  satisfies 


(5.  10) 


B(e2q.v)  =  '''' ®  ” 


where  the  space 


H  vas  defined  in  (3.2)  auid  where 


1 

2 


r(Xj,X2)(v(Xj 


X2,d/2)  +  v(Xj,X2,-d/2))dXjdx2 


(5.11) 


u 


.<1/2 


v(  X  J ,  X2 ,  )  Au(  y ,  q )  dXgdx  jdX2  • 


u  -d/2 

Since  B(e2q,v)  =  0  for  all  veS(2q),  we  can  restrict  v  in  (5.10)  to  be 
of  the  form 


(5. 12) 


d/2 

v(Xj.X2,X2)dx2  =  0, 

-d/2 
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for  a. e.  (x^,X2)€Ci>.  Consequently,  we  arrive  at 

(5.11)'  0  =  R2q(v)  =  I  r(Xj,X2)(v(Xj,X2.cI/2)  +  v(Xj,X2. -d/2)  Icbc^dx, 

*1  f 

*'«  id/2 


2 ) Au( f, q ) dx jdX2dX2 , 


for  all  veS(2q)  where 

(5.13a)  r(Xj,X2)  =  f(Xj,X2)  -  2^^Q^‘^\xyX^,d/2) . 


(5. 13b) 


4u(y.q)  > 

J-0 


Remark  5.3  The  condition  (5.12)  shows,  together  with  Theorems  3.1  and  4.2, 
that  the  bilinear  form  in  (5.10)  can  also  be  understood  as 


B:  X  I - »  R. 

This  implies  in  particular  that  ®2q®^l  qkO. 

Along  the  saune  lines,  B  can  also  be  interpreted  as  a  {C,r)  regular 
form  on  Sj  x  Sg  where 

S[(q):=  S(q)  nH^,  i  =  1,2, 

and  it  is  clear  that 

u(2qj)  -  u(2q2) e  Sj(q2) 


if  0  5qj<q2.  This  is  essential  for  the  adaptive  procedure.  □ 

From  (5.8)  we  have  that 


'0"''2q''l,p 


sup 

V 


IR 


(v) 


.(v) 


2.V 


sup 

V 


,1/2 


(Xi,X2.X3)dXidX2dx3 


8uid  with  (5.11)'  and  Theorem  4.2  we  find  that 
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, 

1.<P  J 


sup 

V 


|r(Xj,X2) |v(Xj,X2.d/2)+v(Xj,X2.  -d/2) 

:d727r~T2 


u 


C[%] 


dx. 


(5.14) 


*  Z 

J-0 


^2^ 


d/2  ,  N 


id/2KJ 


dx„ 


and  the  supremum  is  tsJcen  over  all  O^veH  which  satisfy  (5.12).  It  is 
easy  to  see  that  the  supremum  in  (5.14)  is  attained  on  a  function  v  which 
satisfies  (5.12)  and 


a^v. 


q 


^  J=0 


dv  1 

g— (Xj.X2,±d/2)  =  ? 

w 


tdiere  A2j(x^,X2)  is  completely  determined  by  (5.11)’.  Hence 


(5.15) 

Since 

(5.16) 

we  find  that  in  (5. 14) 

(5.17) 


v(x,,X2.X3)  -Kq.ip]'''*!- 


Xj). 


-d/2 


_  d  1 

3  "  2  (4q+3) 


0  2q  l,fi  q 


12  2 

r  (Xj^,X2)»>  (x^,X2)dXjdX2, 

0 


idiere 
(5. 18) 


.  q  4(4q+3)  ' 

Based  on  the  estimate  (5.17)  we  define  the  error  indicator  function  by 

(5.19)  Tj|  (x  ,x.,)  :*  dC  r^(x  ,  x.,)?»^(x,  .x.,). 

2ql2  q  12  12 
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Remark  5.4.  We  emphasize  that  tj_  is  very  easy  to  compute,  esp>ecially  for 
low  order  models.  In  particular,  we  find  from  (S. 13)  that  for  q  =  0 

(5.20)  tIq(Xj,X2)  =  ^  , 

i.e.  the  indicator  does  not  depend  on  u(0)  and  can  be  computed  a-priori.  □ 
From  (5.14)  and  (5.19)  we  readily  obtain 

and  hence 

(5.21) 


'2q' 


(l> 


'2q"l,<p  Tf, 


e(u(J>,q)). 

0 


Hence  §  based  on  (5.19)  is  an  upper  estimator  and  its  effectivity  index 
6  is  bounded  from  below  by 

Let  us  show  that  8  also  bounds  the  error  from  below.  This  will  allow 
us  then  to  prove  that  8  is  locally  asymptotically  exact  for  auiy  q.  Since 
the  modelling  error  e„  satisfies 


I  I  r(x, 

’^2^ |v(Xi,X2,d/2) +  v(Xj,X2,-d/2) ) 

h) 

-d/2 

q 

<  1  1 

j2Jv(Xi,X2.X3)dX3j 

-d/2  J=0 


for  all  veH  so  that  (5.12)  holds,  we  select  v  =  v<p^  with  v  eis  in  (5.15) 
and  obtain 


(5.22) 


u 

fdv  dv  1 

*  laxj'dxgj 


An  elementary  calculation  shows  with  7^v 


that 


17  vl^Sv^o/lV^vl^  +  6|v|2|7  ®|^}, 

■rV  X  X 

hence 
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(5.23) 


-d/2  P 

"''"L'f  f  K 


dX2dx3. 


-d/2  w 


Id/^  -  2 

I  |3y^l7^v|2  +  6|v|2|7^v|2  +  v"2  1^1  |dXjdX2dx, 


(5.23) 


-d/2  w 


S  3dVDj|r|2_^  »  Q2,rl2  J 


+  dC  irlll^ 
q  0,^ 


where  C  Is  defined  by  (5.16),  Q  by  (4.13)  and 


(5.24) 


and  we  used  the  notation 


D  =  2 


^  (4q+3)  -4 


(5.25) 


'U  ■  f 


|Vxr|^<p^dXidX2,  ir 


i2  [22. 

»0.v  =  J  r  V  dxic 


Now  we  estimate  in  (5.22)  for  c>0 


„  „2  ^  c  „  „2  1  „  „2 

^  2  '’®2,q"l,»>  2e 


If  we  select  >  0  so  that 


(5.26) 


we  get  the  desired  bound 


Using  (5.23)  and  (5.26)  we  find 


(5.27) 


and  we  have 


•  ^ 

n  p  m*  HI  e 

e-  *  ■  1  +  3d‘‘c  D  Q  + - ^ 

0  I  q  ,^,2 


Theorem  5.2.  Assume  that  feT  *  :=  {f|  either  r(q)  *  0  or 

p 

|r*l,  /|r|rt  S6<«o>.  Then  the  estimator  6  baised  on  (5.19)  is  locally 
l.f  0,jp 


(Xj, Kg) -proper  with 


17 


•^1  =  ^0 


f  j  j  j  >  1/2 


(see  (4.10))  and  with 


(5.28)  K_  *  {1  +  3C  D  d^O^  +  Q^)}^^^. 

2  q  q 

Here  and  are  explicitly  given  in  (5.18)  and  (5.24).  If  y  =  1, 

then  Q  =  0  and  the  factor  3  in  (5.28)  can  be  dropped. 


Remark  5.5:  With  the  weight  ^  as  in  (5.5)  and  a  =  1,  we  find  eanily 

that  Q^d^,  0<p<l.  Hence,  for  feT  _  for  some  P  independent  of 

Id  ^ 

d,  the  estimator  S  is  locally  asymptotically  (as  d — »0)  exact.  □ 


Remark  5.6:  From  (5.18)  and  (5.24)  we  see  that  — >1  an  q — for 

f € Tg,  i.e.,  S  based  on  (5.19)  is  also  locally  spectrally  exact.  □ 

Remark  5.7.  The  local  enymptotic  exaictness  of  8  not  only  insures  that  the 

indicator  fimction  ii_  in  (5.17)  gives  a  good  overall  estimate  S,  but  also 

^q 

that  the  local  indicator 


(5.29) 


/  0  0, 
®2q^^*^l*^2^ 


dC 


I 


r  (Xj,X2)dXjdX2 


I I + I I <d^ 


1/2 


is  an  asymptotically  exact  meanure  for  the  local  contribution  to  the  modelling 
error.  This  is  an  essential  ingredient  for  the  adaptive  selection  of  the 
model  order.  □ 


Remark  5.8  Our  analysis  is  obviously  also  valid  for  w  =  (-1,1)  and  we 
present  now  an  example  to  demonstrate  the  shairpness  of  the  bounds  in  Theorem 


5.2. 

For  Q  *  (-1, 1)  X  (-d/2, d/2),  let 

US  select 

f  in  (3.1)  SO  that 

u(x,y)  a  2cos 

[|»)c<,sh  [I 

Then 

q 

q 

«  % 

u(q)  = 

a  2^XjCOS 

J=0  J=0 


T 

The  vector  x  *  (x_, ''^x  )  is  determined  from 

0  q 
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where 


and 


( ^  A  ♦  a  b)x  - 

j 

-ij  “  I  “  I 


a  =  2ir  sinhj— I,  e  =  (1, •••,!) 

.1  .1 


We  find  explicitly  for  the  weight  ^  »  1 

2 

«  -  II  ^ 


'2q' 


‘  «[2cosh[5^]  -  !l\) 


auld  the  estimator 


=  %[*  -  a 


i=l 

Using  a  computer  algebra  system,  we  obtained 


(0(q))2  =  ^/lleoJ?  =  1  +d2—  +  e)(d^) 
2q  1 , p 


where  m  is  listed  in  Table  5.1. 

q 


q 

0 

1 

2 

3 

4 

5 

6 

m 

q 

240 

360 

936 

1768 

2856 

4200 

5800 

Table  5. 1  m^  in  the  eisymptotic  expamsion  of  the  effectivity  index 
Not  only  is  the  estimator  S  asymptotically  exact  as  predicted  in  Thp^i^em 

mm 

5.2,  but  we  observe  that  with  Q  =  0  and  ^  g  ratio  of  (<2^  “  ^ 

2 

(5.28)  sold  n  /m^  in  the  expansion  of  @(q)  is  1,  i.e. 


f  2  2”i 


is  the  best  possible  value  for  p  =  1. 

Let  us  now  analyze  the  error  estimator  for  the  norm  ||e||  ,  =||elL. 

To  this  end  we  consider  the  bilinear  form  (4.15)  and  get 


in 
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(5.30)  Bl'Oa,-'"  -  '<2q‘''> 

With  R2q(v)  given  in  (5.11)'. 

Hence,  using  Theoren  4.2,  we  get 


|R2q(v)| 

Av|| 


Ubar^  the  suprenus  is  taken  over  all  0*v€H2  which  satisfy  (5.12). 

To  estisate  the  supreaus,  we  observe  that  emy  v€H2  can  be  written  in 


the  fora 


(5.31) 


vcxj.xjj.xj)  -  ^]  Y.  *  [a]  Zwk'>‘i-=^> 


k.  ^^l 


trtiere  denote  the  eigenfunctions  ( orthonormal i zed  in  uiu))  of  the 


eigenvalue  problem 


»  0  on  a« 


and  corresponding  sequence  of  ix>sitive  eigenvalues.  Note  that 

RgqCv)  is  as  in  (S.  ID'  and,  as  before,  “  -(4J+1)  r  {XyX^)/d. 

Further,  since  I^q^v)  -  0  for  any  v(Xj,X2,Xg)  *  ^pix^.x^)  =  0  on 

dw,  we  can  oalt  the  second  tern  in  the  expression  (5.31)  for  v  and  find 


1/2 

■*''  *  [a]  Z  ‘’k«'k'’‘i’’^>'=°®(^] 


k,  ^kl 


idiere 


\t  '  v(\*  [t]  ] 


Iisertlng  (5.31)  into  (5.11)',  we  find 


V''’  ■  [a]  f  I  Z  *k«"'V'k‘*i-*2* 

k,eki 


where 
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f 

^  ^(4J+1)  L2q(z)cos(£irz)dz 

J=1  u 

(eukl  a  0  V£).  and  we  estimate 

(R2,(v))^  .  |r,|=3^  J 


idiere 


k^l 
b,..[l-A_J 


til  til  k  d 

»  • 


>dilch  Implies  with  Schwartz’  Inequality  that 


kil  til  k.  til 


Altogether  hence 


(5.32) 


%d)ere 


(R_  (v))*' 

2'’  h!  iav|%  .  L=(„) 

^  L^cn'*) 


E  -  j-  y  V 

q  8-4  Z.  .4 

til  ^ 


Remark  5.9.  For  q  «  0,  we  have  A^^  »  0  and  hence  =  C(4)/(8ii  )  = 

1/720.  For  q>0,  can  be  easily  computed  numerically.  The  bound  (5.32) 
Is  better  than  the  estimate 

^  L  (n  )  w  L  (w) 

tdiich  Is  obtained  from  (5.21)  with  ^  »  1  emd  from 

*®2oll  2  H  ^  !  '®2„l  1  rf  • 

L^(n^)  *  H^di**) 

The  estimate  (5.32)  is  in  fact  optimal,  as  Is  demonstrated  in  Table  5.2. 
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q 

d^llr||%  /  A 

^  2q 

0 

720  +  («d)^/60480  +  O(d^) 

1 

2520  f  (iid)^/221760  +  O(d^) 

2 

10296  +  1287(iid)^/35  +  (?(d^) 

3 

26520  +  9945(ird)^/209  +  Old^) 

2  2 

Table  5,2.  Asymptotics  of  L  -residual  versus  L  -error  for  small  q 
for  the  example  In  Remark  5.8. 


In  all  cases  the  numerical  value  of  the  leading  term  was  found  to  be 

equal  to  that  of  1/E  which  suggests  the  asymptotic  and  spectral  exactness  of 
2 

the  L  a-posteriorl  estimator 


d3/2^ 

q 


2 


□ 


6.  Adaptive  Selection  of  the  local  model  order. 

We  consider  now  the  (?’,q)-model  based  on  an  arbitrary  partition 
f  “  ^*^l^lSlSn  described  in  Section  3,  where  q  =  ^qi»‘‘’*qj^^  vector 

of  model  orders  on  u^.  In  the  previous  section  we  showed  that  the  local 
contribution  to  the  modeling  error  corresponding  to  Wj  cam  be  reliably 
estimated  by  a  calculation  involving  the  local  residuals  only.  This  gives 
rise  to  the  following  simple  feedback  procedure  to  adjust  the  model  orders. 

i)  given  a  parameter  0<X<  1  and  an  order  vector  q,  compute  the 
local  indicators  on  the  domains  according  to  (5.19). 

2 

(6.1)  il)  For  the  nA  of  the  largest  error  indicators  raise  q^  by  one 
and  solve  again. 

iii)  Stop  tdien  Six. 

The  strategy  (6.1)  is  amalogous  to  one  version  of  adaptive  finite  element 
methods.  We  emphasize,  however,  that  (6.1)  is  not  optimal  since  often  q^  heis 
to  be  increased  by  more  than  one  emd  the  computational  solution  of  the  current 
(^,q)-mo<lel  is  costly. 
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Let  us  present  another  adaptive  strategy.  We  start  with  the  observation 
that  for  the  error  e(y, q)  we  have 

(6.2)  B(e,v)  =  R(v)  Vv e H. 

By  Remark  5.3,  we  know  that  e  €  H^,  and,  selecting  <p  suitably,  it  follows 

that  e  depends  only  on  the  local  residual  on  R^.  We  can,  therefore,  obtain 

ain  (aisymptotlcally  exact)  approximation  to  e  on  £1^  :=  u>  x  (-d/2, d/2)  by 

i  d  ““d 

solving  (6.2)  approximately  on  spaces  of  functions  which  vanish  on  £1x5“.  We 

emphasize  that  these  problems  are  completely  decoupled  euid  can  be  solved  very 

2 

cheaply  and  in  parallel.  Therefore,  we  can  assume  that  the  indicators 
Eu*e  known  for  several  values  of  and  all  6 f . 

There  are  many  ways  to  optimize  models  from  the  hierarchy  of  (y.q) 
models.  Any  concept  of  optimality  involves  cost,  e.g.  we  assume  that  the  work 
in  the  numerical  solution  of  the  (?*,q)-model  is  given  by 

n 

(6.3)  W(?>,q)  = 

i=l 

where  the  are  monotonically  increasing  functions  depending,  e.g.,  on  the 
computer  structure,  implementation,  etc.  Then  we  can,  beised  on  a  given  model 
ui?*, q),  find  an  optimal  order  vector  q  by  solving 

n 

(6.4)  min  subject  to  I  nfCqi)  - 

i=l  i 


where  x  is  given  tolerance.  This  problem  has  at  leeist  one  solution. 
Moreover,  using  a  Laingrange  multiplier  A  and  assuming  that  qeIR,  for 
simplicity,  we  find  that  at  the  optimum  necesseu'ily 


(6.5) 


(6.6) 


or,  using  difference  quotients,  that  q  should  satisfy 
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(6.7) 


=  X,  1  S  i  <n. 


Fj(q°PS  -  FiCq!^*"-!) 

This  is  a  weighted  equilibration  condition  characterizing  the  optimal  (f,q)- 

■odel.  Once  we  have  determined  q  from  (6.4)  we  solve  for  uCT*, q)  and 

2 

possibly  repeat  the  process  with  improved  estimates  for  the  The  strategy 

is  capable  of  predicting  directly  the  optimal  distribution  q  and  uses  very 
few  iterations.  Under  suitable  assumptions  the  orders  q^  found  in  this  way 
are  very  close  to  the  best  possible  ones,  hence  the  modelling  is  an  adaptive 
one  (i.e.  a  feedback  procedure  with  certain  optimality  profierties) . 

We  presented  here  only  the  key  ideas  and  refer  to  [41]  for  details  and  a 
more  rigorous  emalysis. 

7.  Adaptive  modeling  in  mechanics 

We  have  shown  in  the  previous  sections  the  main  ideas  of  a-posterori 
error  estimators  and  the  adaptive  modeling.  Although  we  have  addressed  only  a 
roost  simple  model  problem,  ideas  of  this  type  eire  applicable  in  general  (for 
more  see  [ 39 ] ) . 

Various  models  have  different  properties  of  the  solution.  For  example, 
they  have  different  boundary  l2iyer  behavior,  different  singularities  in  the 
neighborhood  of  the  corners,  etc. 

We  emphasize  that  the  estimators  introduced  above  are  not  asymptotical  in 
nature.  They  provide  estimates  for  domains  with  a  thickness  which  is  not 
small.  It  is  essential  that  these  estimates  are  available  for  vairious  norms, 
because  in  practice  the  alms  of  the  modelling  differ. 

Let  us  also  mention  that  the  ideas  explained  briefly  above  are  applicable 
to  many  different  model  -  formulations,  for  exeunple,  to  models  based  on  mixed 
methods  etc.  Various  formulations  also  behave  differently  when  they  are 
descretized  by  the  finite  element  method.  As  example  we  mention  to  the 
locking  problem  [431. 

We  have  assumed  here  that  the  "true  problem"  wets  the  three-dimensional 
problem  for  completely  specified  data.  In  engineering,  however,  often  the 
available  data  have  intrinsic  uncertainties.  These  effects  have  to  be  treated 
by  a  proper  formulation  of  the  true  problem  which  accounts  for  these 
uncertainties,  such  as,  for  example,  a  stochastic  one. 
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